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Abstract 

The scattering of temporally shaped intense laser pulses off electrons is discussed by means of 
manifestly covariant quantum electrodynamics. We employ a framework based on Volkov states 
with a time dependent laser envelope in light-cone coordinates within the Furry picture. An ex- 
pression for the cross section is constructed, which is independent of the considered pulse shape 
and pulse length. A broad distribution of scatted photons with a rich pattern of subpeaks like that 
obtained in Thomson scattering is found. These broad peaks may overlap at sufficiently high laser 
intensity, rendering inappropriate the notion of individual harmonics. The limit of monochromatic 
plane waves as well as the classical limit of Thomson scattering are discussed. As a main result, a 
scaling law is presented connecting the Thomson limit with the general result for arbitrary kine- 
matics. In the overlapping regions of the spectral density, the classical and quantum calculations 
give different results, even in the Thomson limit. Thus, a phase space region is identified where the 
differential photon distribution is strongly modified by quantum effects. 
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I. INTRODUCTION 



The use of chirped-pulse amplification [1] has led to a prodigious advance in available laser 
power. The current records reach several petawatts, and accompanying interest in strong- 
field physics culminates in planned large-scale laser facilities such as the anticipated "Extreme 
Light Infrastructure" (ELI) [2]. The pioneering theoretical studies in strong- field physics 
considered both pair creation in a strong field [3] and the cross channel process, electron 
photon scattering [4-10], dubbed non-linear Compton scattering, where the use of laser 
beams has already been suggested. Since then there has been a wealth of theoretical papers 
and we refer the reader to the reviews [11-16]. In non-linear Compton scattering 

e(p)+e lL (k)^e'{p')+ 1 (k') (1) 

a number of £ photons with momentum k from a high intensity laser, scatter off an elec- 
tron with momentum p. A convenient measure of laser intensity is the dimensionless laser 
amplitude a = eE/mu, with E being the root-mean-square electric field and u the laser 
frequency. The parameter a is a purely classical quantity, representing the work performed 
by the field on the electron in one wavelength. Thus, a is the classical nonlinearity parameter 
[17] and it is related to the ponderomotive potential U p = ma 2 /2. The definition of a can 
be made explicitly Lorentz and gauge invariant [18]. When a becomes of order unity the 
quiver motion of the electron in the laser beam becomes relativistic in a classical picture. 
The spectrum of non-linear Compton scattering has been observed in several experiments 
colliding laser and electron beams, such as low-intensity laser photons (a = 0.01) with low- 
energy (~ 1 keV) electrons from an electron gun [19], a = 2 photons with plasma electrons 
from a gas jet [20] and, more recently, sub-terawatt photons (a = 0.35) from a CO2 laser with 
60 MeV electrons from a linac at the BNL-ATF [21]. Using linearly polarized photons the 
latter two experiments [20, 21] have analyzed the characteristic azimuthal intensity distri- 
butions confirming quadrupole and sextupole patterns for the second and third harmonics, 
respectively. Recently, the energy spectrum of the scattered radiation has been measured in 
an all-optical setup using laser accelerated electrons [22]. 

Probably the best known experiment is SLAC E-144 probing strong- field QED using a 
terawatt Nd:glass laser (a ~ 0.6) in conjunction with high-energy (46.6 GeV) electrons [23]. 
The observation of non-linear Compton scattering has been reported [24] as well as the 
observation of the crossed process of non-linear pair creation, due to the interaction of a 
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Compton scattered high-energy photon with a second laser beam [25]. 

The low-energy limit (in terms of laser frequency) of Compton scattering is Thomson scat- 
tering which is described completely classically [26, 27]. This classical picture is used as the- 
oretical framework for many applications of laser Compton scattering such as X-ray sources 
[28-30] or diagnostic tools [31]. A convenient parameter to distinguish the two regimes is 
the quantity 

s e - m 2 2£k ■ p 

ye = — — = — a - , 2 

where s e = (p+£k) 2 expresses the center of mass energy squared for the generation of the £th 
harmonic in a Lorentz invariant manner 1 . The £+ 1 — >■ 2 process is kinematically equivalent 
to the scattering of one photon with momentum £k off an electron with momentum p, thus it 
appears as a (pseudo) 2 — > 2 process. The Thomson regime is recovered for yt <^ 1, while for 
Hi > 1 one finds striking differences to the Thomson scattering. The electron recoil during 
the scattering may be quantified by the Lorentz-invariant quantity t — {p — p') 2 which is in 
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the range < — t < rn 2 j^, i.e. in the Thomson regime —t/m 2 <C 1 holds. 
A quantity measuring non-linear quantum effects, 



ey^Fggg 1 

XR = § = o^i, 3 

m 6 2 



has been introduced in [4, 8]. It measures the work done by the field over the Compton wave- 
length mr l in the rest frame of the initial electron where the four-vector p^ = (m, 0,0,0). 
Introducing the critical field strength [32] Eg = m 2 /e = 1.3 x 10 18 V/m, this may also be 
written as xr — E+/Es, where E± is the rms electric field strength in the electron's rest 
frame. The parameter xr combines nonlinearity and quantum effects, xr is °f the order 
of unity if both a and y 1 are of the order of unity. Thus, the corrections to the classical 
description (Thomson scattering) are important if either (i) an ultraintense high-energy pho- 
ton pulse, e.g. produced by an X-ray free electron laser, interacts with low-energy electrons, 
or (ii) a multi-GeV electron beam is brought to collision with an optical high-intensity laser. 
The latter scenario is similar to the SLAC E-144 experiment but with a higher value of a. For 
the 50 GeV SLAC beam in conjunction with a counterpropagating optical laser {uj ~ 1 eV) 
one has yt ~ 1. These parameters will be used mainly below for numerical calculations. The 



1 p a,n k are four- vectors, thus k ■ p denotes a scalar product; we employ units with h = c = 1. 
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FACET project [33] at SLAC envisages investigations within such kinematics in line with 
(ii). 

Ultraintense lasers use short pulses (few fs, few laser cycles) requiring a proper treatment of 
the laser pulse structure. Rich substructures of the scattered photon spectra were predicted 
within the classical picture [34-37] of Thomson scattering. These substructures have not 
yet been confirmed experimentally. The effect of radiation back reaction on the spectra 
was studied in [38] and found to be important. Only a few publications address quantum 
calculations in pulsed fields for scalar particles [39] and for spinor particles [40, 41]. In [42], a 
connection between the emitted angular spectrum in non-linear Compton scattering and the 
carrier envelope phase in few-cycle laser pulses was established. In a related field, electron 
wave-packet dynamics in strong laser fields has been studied in e.g. [43, 44]. 
In our paper we calculate the emitted photon spectrum in non-linear Compton scattering 
using a generalized Volkov solution with temporal shape in light-cone coordinates. For this 
purpose we first focus on the structure of the Volkov wavefunction in a pulsed laser field. 
The aim of our study is to compare the QED calculations for non-linear Compton scattering 
with results from a classical calculation, i.e. Thomson scattering. 

Our paper is organized as follows. In section II we present Volkov states in pulsed laser 
fields. Section III continues with the calculation of the matrix element and the transition 
probability. A slowly varying envelope approximation is discussed. Numerical results are 
presented in section IV. We discuss various limiting cases of our general results, including 
monochromatic plane waves and the Thomson limit. As a main result, a scaling law is pre- 
sented, connecting the Thomson spectrum with the Compton spectrum. In the Appendix we 
summarize the kinematics in light cone coordinates and present the Fourier transformation 
of the Volkov state. 

II. TEMPORALLY SHAPED VOLKOV STATES 

A strong laser field may be considered as a coherent state of photons \C), characterized by 
the polarization and momentum distribution C^{k), if the depletion of the laser photons 
from \C) by an interaction process with electrons is negligible, i.e. for any relevant scattering 
process S = (out;C'\S\in;C) with C — C is valid, where in and out are particle number 
states without coherent parts [45]. Then, it is possible to work within the Furry picture [46], 

4 



where the interaction of an electron with the classical background field A^(x), which is the 
Fourier transform of C M (fc), is treated nonperturbatively and solutions of the Dirac equation 

{%$ -e£- m)il>(x) = (4) 

are utilized as basic in and out states for the perturbative expansion of the S matrix. For 
background fields in the form of plane waves, closed solutions of (4) can be found, 

where the free Dirac spinor for momentum p and spin s fulfills Qzi — m)u PtS = and is 
normalized to u p ^u p ^ s = 2m5 ss '. The phase is the classical Hamilton Jacobi action 

S p (x) = -p-x + f(k-x), (6) 

k-x k-x 

with / = fx + f 2 , fi = — J d0^^ and f 2 =f Equation (5) represents the famous 

Volkov states, whose perturbative expansion in terms of interactions with the laser field 
is depicted in figure 1 of [45] (The expansion parameter, i.e. the coupling strength at the 
vertices, is a defined below). 

For the vector potential we use a real transverse plane wave 

= A g{k ■ x) (e^ cos£ cos k ■ x + e 2 sin£ sin k ■ x), (7) 

modified by an envelope function g and fulfilling A ■ k — and k ■ k = 0. The parameter £ 
determines the polarization of the laser: It is linearly x (y) polarized for £ = (£ = 7r/2) and 
circularly polarized for £ = ±7r/4. For other values of £, the laser is elliptically polarized 
[47]. The vector potential is normalized such that the mean energy density or the energy 
flux (E 2 ) oc — A^A^ = g 2 Al/2, where (. . .} means averaging over the fast oscillations of 
the carrier wave, is independent of £, but the dimensionless laser amplitude a, as defined 
in the introduction, is time dependent. A time independent laser strength parameter may 
be defined by the normalized peak value of the vector potential ao = eAo/m (with this 
definition, a 2 = a%/2 for g = 1). The vector potential (7) can also be cast into a complex 
form 

A» = Ao9{k ' X) {e%e' ik - x + e"_e +ik - x ) = ^-B^ (7') 
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with the complex polarization vectors e± = cos£e^ ± ism^ with e + ■ e_ = —1, e± • e± = 
sin 2 £ - cos 2 £ and the definition B» = e%e~ ik - x + ete +ik - x . 

In what follows, the temporal pulse shape will often be chosen as a cosh pulse 

g{k-x) = (8) 

cosh I — 

with width <r, or a Gaussian 

^•x)=exp(-^-^. (9) 



2a 2 

Besides these special cases, any other smooth function g which depends solely on k ■ x is 
possible. 

Since the vector potential A 11 depends only on k ■ x, it is convenient to work in light-cone 
coordinates with k-x = ux + (see Appendix A). In these coordinates, the Volkov wavefunction 
(5) reads 

ip p , s (x) = C p (x ± ,x-,x+) ^p= , (10) 

y/ lbj v 

C p (x ± ,x^,x + ) =\l + d p g{x + ) (^_e i0JX + + ^ + e~ iujx +) 1 e -i(p+*-^-*+)-Hp±-*L+i/(*+) (n) 

with dp = aom/(4:k ■ p). For many purposes it is sufficient to consider only the function C p 
since it contains the relevant information on the interaction of the electron with the laser 
pulse. The real part 2 of the scalar projection S —^trC p (x) = exp{iS p (x)} is visualized in 
Fig. 1 in the frame where the electron is initially at rest. The scalar projection is essentially 
equivalent to probing the state tp PtS with u PyS and an average over the spins s 

In that frame the (free) electron wavefunction outside the laser pulse behaves as oc e~ ip ' x = 
e -imt rpj^ e g ec ^ £ laser pulse is a local deformation of the electron wavefronts due 
to the build-up of an effective, time dependent momentum g M (x + ) = p^ + g 2 (x + )b p k fl with 
b p = m 2 a^/Ak ■ p. The momentum q^(x + ) achieves its maximum at the center of the laser 
pulse x + = 0, depicted as diagonal straight line in Fig. 1, where it coincides with the usual 
quasi momentum in monochromatic plane waves. Thus, inside the laser pulse, especially 



The imaginary part gives no further information; it has a shifted phase as compared to the real part. 
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Figure 1. Contour plot of ^txC p (x) in position space in the z — t plane. The laser pulse with 
cosh profile is located between the two dotted lines. Left (right) top panel: Circularly (linearly) 
polarized laser pulse with ao = 1.5 and a = 20, Bottom panels: circular polarization for ao = 0.5 
(left) and ao = 3.0 (right). 

for x + = 0, the fully dressed electron wavefunction behaves as oc e~ iq ' x = e - % ( m + b P u >*) t +' lb P u * z ^ 
i.e. the electron wavelength changes and the wavefronts become tilted. Both effects are 
proportional to the ponderomotive potential, i.e. oc b p u* = ma^/A = U p /2 , where u± = u- k 
is the laser frequency in the initial electron rest frame and u = p/m. In Fig. 1, = 300 keV 
is chosen which can be achieved, for instance, with laser photons of 1.5 eV energy colliding 
head-on with 50 GeV electrons, or a 15 keV X-ray laser beam with 5 MeV electrons. The 
electron wavefunction changes its behavior from the free case to the fully dressed case over 
iV = am j (271U+) oscillations of the free electron wavefunction, i.e. iV = 5.5 for the parameters 
employed in Fig. 1. For lower values of u;*, the behavior of the electron wavefunction changes 
slowly over many oscillations, e.g. for = 200 eV (35 MeV electrons colliding with 1.5 eV 
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Figure 2. Contour plot of the spin flip contribution to the Volkov wave function by the tensor 
projection T 13 — iT 23 . Same parameters as in Fig. 1. 

photons) and o = 20 one finds iV = 8100. For a linearly polarized laser, additional ripples 
appear in the interaction region due to the oscillating terms in the phase. For non-head-on 
geometries, similar ripples are also present for circular polarization. The smoothness of the 
pulse envelope g ensures the smoothness of the Volkov wavefunction in the transition from 
the field-free regions to the laser pulse. 

Vector projections with an odd number of Dirac matrices, such as V M =|tr •y^C p (x) or 
A? =|tr 7 5 7 M C p (x), vanish identically and are, therefore, not useful for characterizing ip PtS . 
The pseudoscalar V = |tr7 5 C p (x) also vanishes. The antisymmetric tensor projections 
jw = |tr a^C p (x) = id p (k u A» - WA U ) exp{iS p (x)}, where o» v = is the spin ten- 

sor, allow for a further characterization by a 01 , a 02 , a 13 , a 23 for circular laser polarization and 
for for linear polarization. These tensor projections are nonzero but only inside the 

laser pulse. They mix contributions with different spin orientation and are therefore propor- 
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Figure 3. Feynman diagram for non-linear Compton scattering as the decay of a laser dressed 
Volkov electron state. 

tional to a combination of Up^ipp^ix) and u v ^v,\ ( x )i i- e - the spin-up wavefunction contains 
contributions with spin-down and vice versa. From the structure of the Pauli interaction 
term a^F^, where = d^A u — d u A^ is the electromagnetic field strength tensor, one can 
infer that a 01 (a 02 ) corresponds to the interaction of the electron with the x (y) component 
of the electric field, and a 13 (c 23 ) corresponds to the y (x) component of the magnetic field. 
From this correspondence and by inspecting (7) it is easy to understand why some projec- 
tions are zero for linear polarization while they are nonzero for circular polarization. As an 
example, the tensor projection T 13 — iT 23 is shown in Fig. 2. 



III. CALCULATION OF THE MATRIX ELEMENT 
A. The S matrix 

The interaction of the Volkov electron ey with photon modes different from the laser field are 
treated by perturbative S matrix expansion. The Born approximation of the matrix element 
for the emission of one photon, i.e. non-linear Compton scattering ey(p) — > ey(p') + j(k'), 
is depicted in Fig. 3. Using Feynman rules [48], the S matrix element for such a process is 
given by 



S ft = (p',s';k',e' x ,\S[A}\p,s} = 
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which reads in light-cone coordinates, suppressing spin (s, s') and polarization (A') indices 
from now on, 

S fi = N J d 4 xu p/ (l + d p ,$%)f.'{l + d p 0)u p e i{Sp - s p' +ik '- x) (14) 
at r 

= Y d 2 x±dx + dx_T(x + )e iH ( x +> x -' x ^ (15) 
with iVo = -ie/y/2io'2E p 2E p/ and 

T(x+) = + g e tuJX+ + g e~ iu)X+ &\ + g 2 ST 2 + g 2 e 2iux+ ST 2 2 + g 2 e- 2twx+ 2T 2 2 , (16) 

where 

= u p ,f'u p , (17) 

^ = v (^'^^' + «p, (18) 

<5q 2 = 4(/c • e')dpdp'U p iftu p , (19) 
^ = d v d v ,u v , {i^i'^Q Up. (20) 

Due to ^ ftffy = — 2(e' ■ fc)(e T • e T )^, one finds ^ 2 = for circular polarization. Further- 
more, 

H(x+, x_, x_l) = — £y + ik' ■ x 

= (p' + k'-p)-x + f(x + )-f'(x + ) (21) 

with /' = f(p -+ p') = hip 1 ) + f 2 (p') and 

k-x 

fi(x+;p) = / d05-(0) [p-eicos^ cos0 + p -easing sin0], (22) 

"-■pi 

00 

TYl 2 CL 2 f 

f 2 (x+;p) = -y^J d0/(0)[cos 2 £cos 2 + sin 2 £sin 2 0]. (23) 

Inspecting Eq. (21), it is obvious that the dependence of H on x_ and Xi is trivial and the 
integrations over these variables in Eq. (15) can be done analytically. As a result, energy- 
momentum conservation is imposed on the components P + and Pj_, and the exponent 

H + {x + ) = l -{k'_ + p'_-p,) x+ + f( x+ ) - f( x+ ) (24) 

remains. Due to the non-trivial pulse dependent structure of H + , the x + integration does 
not yield another conservation law. Thus, the frequency of scattered photons uj' is not fixed 
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by energy and momentum conservation as a function of scattering angle 9 and remains as 
independent parameter. Including the x + dependence of F(x + ), some rather complicated 
functions of oj' emerge 

oo 

**n = J dx + g M (x + )expi{H + (x + ) + Nux + }. (25) 

— oo 

With these definitions, the S matrix element can be written as 

S fl = (2n) 3 5 2 (k' ± + p' ± - p ± )S(k' + + p' + - P+ )N JZ (26) 

with 

JK = ^V ° + ST^l + &\s*\ + ^V 2 + ^2^2 + ^-W-V (27) 

The integrals srf™ are numerically convergent for M > 1 due to the presence of the pulse 
function in the integrand, rendering the range of integration practically finite. The integral 
£*o°, however, contains a divergent part and must be regularized. A possible method has been 
proposed in [41], where one multiplies the integrand with a convergence factor e~ £][X+ \ e > 0, 
and performs an integration by parts. The result is 

oo 

< = -^- J dx + d( ^~ fn> e W {iH + (x + )} + 4jW)-rm hm , (28) 

— oo 

with P_ = k'_ + p'_ — p-. In (28), the first part is now convergent and the second part 
is proportional to a 5 distribution with support at uj' = 0. The latter contribution can be 
neglected in our analysis for uj' > 0. This regularized version of srf® will be used in the 
subsequent numeric calculations. 



B. Slowly Varying Envelope Approximation 

The calculations are simplified upon utilizing the slowly varying envelope approximation 
(SVEA) of the phase of the srfff functions. This approximation scheme is suitable for long 
pulses with a ^> 1. Typically a is proportional to the number of laser oscillations under 
the envelope. For f\, which is proportional to g (see Eq. (22)), an integration by parts is 
performed, yielding 

j d<fig((f)) sin0 = — g(<f>) cos0 + j d0-^cos0, (29) 
j d(f>g((f)) cos = g{4>) sin — j d(f)-^ sin <p. (30) 
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SVEA basically means neglecting the second terms containing the derivative of the pulse 
shape dg/ d0 because it is 0(1/ a) smaller than the first term. For f 2 (oc g 2 , cf. Eq. (23)) 
we use 



/ 
/ 



d0g 2 (0) cos 2 ■ 
d0g 2 (0) sin 2 i 



&<pg 2 (4>) + ^g 2 ((p) sin cos 0, 
#^ 2 (0) - ^9 2 (<P) sin cos 0, 



(31) 
(32) 



which becomes particularly handy if f d(pg 2 is known analytically, such as for the sech 
pulse (8), where J d0cosh~ 2 0/cr = crtanh(0/(j) + const, or the Gaussian pulse (9), 
J d0exp(— 2 /2<t 2 ) 2 = v / 7T(Terf(0/(T)/2 + const, where erf(x) is the normalized error func- 
tion. Finally, the SVEA result for the phase reads 



ma Q 



fi = g(x + ) [p ■ e\ cos £ sin oox + — p ■ e 2 sin £ cos oox + ] , 



k ■ p 



f2 



m 2 a\ 
4k ■ p 



UIX + 

J d0g 2 (0) + g 2 (x + ) coswa; + sincjx + (cos 2 £ — sin 2 £) 



(33) 
(34) 



generalizing the approximation scheme of [40] to linear laser polarization. 
Even for short pulses, such as for a = 5 meaning that there are about 5 laser oscillations in 
the pulse, i.e. the pulse length is ps 15 fs for A = 800 nm, SVEA is quite a good approximation, 
see Fig. 4 for selected examples. 



C. The spectral distribution of scattered photons and the cross section 

In the standard formalism, scattering experiments are thought of as constant streams of 

particles interacting. Consequently, the square of the S matrix contains a factor T which 

originates from the square of the energy- momentum conservation which is interpreted as 

8(Pi — Pf) 2 — > j^yi5(Pi — Pf), with the volume V and interaction time T which are both put 

to infinity. On the purpose of rendering this quantity finite, usually the differential rate per 

\Sf l 2 

unit time and unit volume dw^f = y/, dH is considered, where did denotes the final state 
phase space. Here, however, the interaction is happening only within a finite time interval. 
Because of lacking one 8 distribution, the square of the S matrix now reads 

\S fl \ 2 = (27i) 3 V5(p' ± + k' ± - p ± )5(p' + + k' + - P+ )\N ^\ 2 , (35) 

where the dependence on the interaction time is contained in ^# and is finite. Thus, it is 
not necessary to define a differential rate per unit time. An appropriate observable is the 

12 



circular polarization 



linear polarization 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

uj'/co 



Figure 4. Comparison of the SVEA (black dashed curves) and the full numerical results (solid 
curves for 8 = (red), 9 = I/7 (blue) and 9 = 2/7 (green)) for the real parts of the functions £/q 
and for a = 5. Parameters are ao = 1.5, U) = 1.5 eV, 7 = 10 5 . Left (right) panels are for 
circular (linear) polarization. 



emission probability of photons per unit volume and laser pulse 



diV 



I Q |2 

V 



dn, 



(36) 



which has as classical analog the spectral density of scattered photons in Thomson scattering 
(cf. e.g. [26, 37]) 



d 2 N, 



classical 



IjJ 



16tt 3 ~ 
< I dr^(r)e ifc '- x(r) , 



(37) 
(38) 



where u^(r), x M (r) are the classical velocity and orbit from a solution of the Lorentz force 
equation for a spinless pointlike charge, and j^ik') is the retarded Fourier transform of the 
electron current. The notion of Thomson scattering is specified to mean this particular 
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calculation scheme. A quantum spectral density is given by the Lorentz invariant expression 

d 2 N ss , x , _ eV 1 2 

dTdfi " Te^i^v 1 ^' 1 ' (39) 

which depends on spin and polarization indices. Averaging over the spin of the incoming 
electron and summing over the spin of the outgoing electron and the polarization of the 
outgoing photon yields a quantity which is directly comparable to the classical spectral 
density 

d 2 iVq Uantum _ 1 ^-v v d 2 N ss t\i 

s,s'=l A' 



du'dQ 2 ^ ^ du/dfi 

s,s'=l A' 

Now we construct an invariant cross section by dividing Eq. (40) by the normalized number 
of photons Nl in the laser pulse, i.e. 

d 2 a 1 d 2 A" quantum 



du'dQ N L du'dQ ^ 
with Nl = JZo dt( ?' where S = E 

x B is the Poynting vector of the laser field, derived 

from the vector potential (7), yielding 

oo 
— oo 

with dtcosh _2 0/cr = 2a/u and dtexp (— (f 2 /2a 2 ) 2 = ^Txajuj for the pulse shapes 
(8) and (9), respectively. Using this definition, the total cross section is independent of the 
pulse shape function g and the pulse length a. (A different definition of the cross section 
without this property has been proposed in [41].) This has been checked numerically for 
a — > by a comparison of = J du' d f,^ n with the differential Klein-Nishina cross section 
[48], or (Ttot with the total Klein-Nishina cross section. In particular, in the limit y 1 — > we 
obtain the total Thomson cross section Ut = 665.25 mb accurately, as exhibited in Fig. 5 
for three different pulse shapes and pulse lengths. 

IV. DISCUSSION 

A. Monochromatic Limit 



In the famous case of monochromatic Compton scattering, the frequency of the scattered 
photon is uniquely defined by the scattering angle. For a finite temporal laser pulse, however, 
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Figure 5. Total cross section for Compton scattering normalized to the Thomson cross section cry. 
Red curve: Klein-Nishina cross section. Symbols: numerically calculated cross section in a pulsed 
laser field with ao = 0.001. Blue stars: a = 20 for a cosh pulse; green circles: a = 30 for a Gaussian 
pulse; black diamonds: a = 100 for a Gaussian pulse. 

this tight relation is lost. As outlined in subsection III C, there is a distribution of the emitted 
photons for a fixed angle. As an example, we exhibit in the left top panel of Fig. 6, the 
spectral density d 2 N/ du' dQ as a function of w' = w'/^i classical f° r nxe d ^- The vertical 
thin lines depict the positions of the harmonics for a monochromatic plane wave with infinite 
duration and the same value of ao, given by [48] 

£k-q 



(43) 



.quantum ( ? + ^) • n' ' 

introducing the intensity dependent quasi-momentum of the electron q^} = p^} + b p (,)k_ 
with b p (,) = m 2 a1/Ak -pW and the dressed mass-shell relation q 2 = q' 2 = ml = m 2 (l + Oq/2). 
Note that only p^}, the conjugate momentum to x + , is modified by an intensity dependent 
contribution, i.e. q+ = p^} and = p^. The integer £ labels the individual harmonics, 
which are not equidistant in general. 

In a pulsed laser field, each harmonic consists of a bunch of spectral "lines" (or subpeaks) 
visible in the top panels of Fig. 6 with a certain width Au' e determined by the minimum and 
maximum values of intensity in the laser pulse. The high-energy tail of each harmonic bunch 
is given by u' e (ao — > 0), and is produced at the edges of the laser pulse. The low-energy edge 
is given by oo' e (a ) and accounts for the maximum red-shift at the center of the pulse. Thus, 
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the spectral width of each harmonic t is given by 

Au' e = u' e (a -> 0) — io' e (a ) 

b k-n' (44) 
= u' e {aoM{ao -> °) jj fc • 

The number of subpeaks in a bunch is proportional to the pulse length a and the intensity 
Oq. The highest subpeak takes its maximum value at a higher frequency oj' than predicted by 
(43), thus, at a smaller intensity-dependent red-shift than the monochromatic harmonics due 
to a lower average do- Hence, one could say that this maximum is blue-shifted as compared 
to the monochromatic plane wave. 

Increasing a from 20 to 50 does not lead to an accumulation of spectral weight at the non- 
linear Compton frequencies as could be expected naively. The number of subpeaks increases 
but the average shape of the harmonic bunch is more or less the same for a = 20 and 50 with 
the same spectral width. In fact, to obtain the monochromatic limit, it is not efficient to 
take simply the limit a — > oo. A method with better convergence is to introduce a flat-top 
area in the pulse. This however, introduces a second pulse length parameter: The total 
pulse length now consists of the rise "time" a and the flat-top "time" r. The flat-top part of 
the pulse is parametrized as <7fl a t(0) = ©(0 + ht)Q{iit — 0), where a factor n is introduced 
so that r is comparable to the Gaussian and cosh widths a in terms of laser oscillations 
under the envelope, and 6(0) is the Heaviside step function. Then, the complete pulse is 
parametrized as 

# to t(0; <J, t) = 0(0 + itt)Q(itt - 0) + g((f) - 7tt)0(0 - ttt) + g((f) + ttt)0(-0 - ttt).(45) 

The spectrum converges rather fast to sharp peaks centered at the non-linear Compton 
frequencies upon increasing r from to 30 while keeping a = 20 constant, as seen in the 
bottom panel of Fig. 6: The strengths are located at the sharp non-linear Compton energies. 
The remaining wiggles around the non-linear Compton energies vanish upon increasing r 
further. 

In the monochromatic limit r — > oo, the rising and trailing edges of the pulse shape function 
become unimportant, i.e. g — > 1, and the function H + in (24) reduces to 



H + = \ ( k'_ + q'_ — q_ )x + + a,\ sin tux + — a 2 cos ujx. 



(cos 2 £ - sin 2 £ ) sin 2ux^ 



(46) 



with ctj = mao(ej • p/k • p — • p'/k ■ p') and re-identifying the electron quasi-momenta 
q^} . Upon plugging (46) into (15) and expanding into a Fourier series, one obtains a fourth 
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Figure 6. Top panels: Spectral density d 2 N/ dui 1 dQ as a function of the scaled frequency w' = 
u ' /^i classical f° r a cos ^ pulse with cr = 20 (left) and 50 (right). Bottom panel: Spectral density for 
a flat-top pulse with cosh edges, r = 30 and a = 20. In all panels ao = 1.0, 7 = 10 5 , u) = 1.5 eV, 
0=1/7 and </? = 0. The thin vertical lines depict the non-linear Compton energies defined in 
Eq. (43). 

energy-momentum conservation by integrating over x + , yielding q_ + £k_ = q'_ + k'_ . The 
four energy-momentum constraints together lead again to Eq. (43). 

The individual harmonics, consisting of a multitude of subpeaks, begin to overlap if the 
lower edge of the £ + 1st harmonic coincides with the upper edge of the £th harmonic, i.e. 



u' e (a -)> 0) > Lu' £+1 (a ). 



(47) 



This happens always for sufficiently large values of ao and £. The notion of individual 
harmonics becomes inappropriate, as one rather observes a continuous spectral distribution. 
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Figure 7. Different recoil parameters y(oj') and yi for I = 1, 2, 3 as a function of the scaled frequency 
w' = w'M.classicar 

B. Comparison with Thomson scattering 

There are different bookkeeping parameters for the characterization of the Thomson regime 
as limiting case of the presently considered scenario. One parameter is yi introduced in 
Eq. (2). An alternative would be to employ the outgoing momenta instead of the incoming 
ones, defining y = (s — m 2 )/m 2 with s = (p' + k') 2 . When four-momentum conservation 
holds, both definitions coincide (since k' and p' both depend on £) and s coincides with the 
usual Mandelstam variable s. However, this is not the case here. These recoil parameters 
are compared in Fig. 7. The parameter y is a function of oj', as it depends on uj' through k' 



p'-k' 2 (p-k'){p-k) 2 7 V(1 



n 



(3) 



77T 



m 2 (p ■ k — k' ■ k) 



(48) 



which diverges at cu c 



1— nn' 



, defining the boundary of phase space. Thus the physical 



phase space is given by < 9 < tt, < (f < 2tt and < oj' < oj'^. An interpretation of the 
phase space boundary will be given in subsection IV C. 

To relate the Compton amplitude with the classical Thomson counterpart it is instructive 
to consider the phase exponential, e.g. H + ± ojx + in &/±i, cf. Eq. (25). For the sake of 
simplicity, a backscattering head-on geometry with a circularly polarized laser is assumed 
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in this subsection. Then, after using some light-cone algebra, the phase reads 



H + ± cux + = 



(k'_ +p'_ -p_)±u 
k' ■ p 



x 



pi 



±00 



x + + 



+ 



k'-k 



2„2 



ma, 



o 



2k -p 2k 



m 2 a 2 \ I" 
2h-p>) J 



d0 O 2 (0) (49) 



2„2 



ma. 



o 



ck/>/(0). 



(50) 



(k-p)(k-p') 2 

Momentum conservation implies n + ■ p' = n + ■ p — n + ■ k' (see (A7)). For n + • F < n + ■ p, 
the leading term 

k! ■ p 



x + ±0J 



n + -p 



+ 



k! ■ k m 2 a\ 



d(pg (0) = k' ■ x(t) ± ujx_ 



(51) 



{k-pf 2 

agrees with the corresponding expression obtained in a classical calculation for Thomson 
scattering (cf. [37], see also Eq. (38)). The frequency of back-scattered photons in monochro- 
matic plane waves is obtained from (43) by neglecting ik w.r.t. q in the denominator, i.e. 

k ' q (52) 



U) 



^.classical 



q ■ n 



In Thomson scattering, the harmonics are always equidistant. A series of plots showing 
the transition from Thomson to Compton scattering is exhibited in Fig. 8. The deviations 
between Thomson and Compton scattering are: {%) a non-linear red-shift in frequency and 
(ii) a slight modification in the amplitude starting notably at y\ = 0.12. It is obvious that 
the red-shift is much more pronounced at higher frequencies. Figure 8 quantifies the well 
known fact [48] that Compton scattering turns into Thomson scattering in the low-energy 
limit. For the chosen parameters (ao = 1.0, oo = 1.5 eV) the differences become significant 
for 7 > 10 4 . Very drastic differences are obvious for 7 = 10 5 (bottom right panel of Fig. 8). 



C. Scaling properties of the spectral density 

The classical and quantum spectral densities for arbitrary pulse shapes are connected by the 
scaling law 

d -^classical/ / n\ ^ -^quantum/ / n\ /ro\ 

-d^dn-t"'^ du/dQ {x ^ e) (53) 

with the two scaling factors r\ and x which are determined by the monochromatic results. 
The frequency scaling factor \ is given by 

2 

= ^.quantum = Tl' ■ U + Tl' ■ = q ■ n' 

^classical n> ■ U + U> ■ 71 ( ^_ + £ « ) (q + M) ■ U>' 

\4n-u m J 
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Figure 8. The photon spectrum as a function of the scaled frequency w' = w'/w^ classical (0) for 
7 = 10 2 ,10 3 ,10 4 ,3 x 10 4 ,6 x 10 4 ,10 5 (i.e. yi = 0.0012,0.012,0.12,0.35,0.7,1.2) from top left to 
bottom right for 9 = I/27 and (p = for a cosh pulse shape. The upper (lower) blue (red) curves 
are for a quantum (classical) calculation of Compton (Thomson) scattering. The vertical gray lines 
mark the positions of the monochromatic harmonics Eqs. (43) and (52), respectively. The other 
parameters are qq = 1.0, ui = 1.5 eV, a = 20 and £ = 0, i.e. linear laser polarization. 
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A continuous effective £ cS has to be used, which follows from the inversion of w^ quantum , 
yielding 

, , ,s g ("' ■ * + ■ ^ Q-k' ,,,, 

^ 1 = „■«-„-■,.* = o-t-if-f (55) 

m 

which simplifies to % = 1 — k' ■ k/p ■ k = 1 — k' + /p + . The scaling of the frequency naturally 
also includes the scaling behavior of the phase space factor which is proportional to u' 2 . 
The scaling factor 77 describes the scaling of the differential probabilities defined by 

(1—1 d(Xq Uan t um \ / f _2 d(T c j ass j ca [ \ fKC\ 

^quantum I I ^classical J i V 3 ®) 

where the differential cross sections do" quantU m and do" classica i are the monochromatic plane- 
wave cross sections, yielding for circular polarization [37] 

X " (57) 

a 



where x = (1 — x)/Xi V* — ■ p/m 2 and £ c ff has to be used instead of I everywhere. The 

l+x' 



other definitions are & t = J 2 +1 (z) + J}_ x (z) -2J 2 (z), &t = -8Jf(z)/a 2 +2£ e , Zt = &e+j^£t 



and z = 2£y \J yl^ ~ yl ) ' are Vessel functions of the first kind. In the limit a — > 

one gets 

x 2 1 

lim = H — — (58) 

ao^o ' i + x 2-4^(l- ^) v ; 

which is a good approximation for a < 1. For linear laser polarization we expect a similar 
relation to hold, but with the appropriate linearly polarized monochromatic plane-wave cross 
section instead. 

The scaling function x is related to the momentum transfer from the incoming electron to the 
outgoing electron via p' + = XP+- Thus, x is the fraction of p + momentum transferred from 
the incoming electron to the outgoing electron p' + . Furthermore, the fraction of momentum 
transferred to the photon is k' + — (1 — x)P+ thus, 1 — % is another measure of the electron 
recoil, x i s a monotonic decreasing function of 10' . The point x{ u ') = corresponds to 
p' + = and k' + — p+, i.e. the total amount of momentum is transferred from the electron to 
the photon. Further increasing u' would render x as well as p' + negative. This defines the 
boundary of the physical phase space for the outgoing particles. When p' + < 0, then, due 
to the free particle dispersion relation (A6), also p'_ < 0. This, however, would lead to a 
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Figure 9. Left panel: Different components of the outgoing electron momentum p' x as a function of 
frequency uJ for backscattering geometry in the rest frame of the incoming electron. Shown are p'_ 
(red, dash-dotted), p\ (green, dashed) and p' Q = E p i (blue, solid). The physical phase space has its 
support at x > 0, i.e. for < u'/m < 0.5 = co'^/m in this case. Right panel: Maximum frequency 
cj'oa as function of scattering angle 9, where 9 = denotes the backscattering direction. 

negative energy E p i due to E p / = (p' + + p'_)/2 < 0. Consequently, there exists a maximum 
frequency u'^, defined by xi^'oo) = ; 

, p ■ k 7717(1 — n ■ (3) 



n' ■ k 1 — n • n' 



(59) 



which can also be obtained as the limit lim u' e = u'^, with u' e from Eq. (43). This also 
coincides with the singularity in y, see Eq. (48). For the backscattering head-on geometry 
we obtain u'^ = m/2 for electrons initially at rest and u'^ = 7717(1 + 0)/2 ~ E p for 
ultrarelativistic particles, i.e. the maximum backscattered frequency is determined by the 
energy of the incoming electron. The momentum components p' + , p'_ and E p > are depicted 
in the left panel of Fig. 9 as a function of u' in the electron rest frame. The right panel of 
Fig. 9 shows the dependence of u'^ on the scattering angle 6. It takes its minimum at the 
backscattering direction 6 = and goes to infinity in the limit 9 — > 7r, i.e. forward scattering. 



D. Further differences between the classical and QED calculations 

As demonstrated in the preceding subsections and quantified by the scaling law, the main 
difference between the spectral densities of Thomson and Compton scattering is the proper 
treatment of the electron recoil in the latter one. Additionally, there is another regime where 
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Figure 10. Comparison of the quantum (upper) and classical (lower) spectral distributions 
d 2 N/ doj' dQ as a function of oj' for fixed angles 9 = I/7, tp = 0. The overlap regions of dif- 
ferent "harmonics" are highlighted as gray shaded areas. Parameters are a = 50, uj = 1.5 eV, 
a = 2.0, 7 = 80, i.e. y < 2.5 x 10~ 3 . 



the quantum description goes beyond a classical calculation, even if y <C 1, where the total 
Thomson and Compton cross sections are equal in leading order. This happens in regions of 
phase space where individual harmonics are overlapping (see Eq. (47)). There, the subpeaks 
in the quantum calculation show completely different patterns in comparison to a classical 
calculation, see Fig. 10. Consequently, the scaling law may not be applied where harmonics 
are overlapping. For a better orientation, the spectral ranges of the individual harmonics 
are marked in Fig. 10, where the lower (upper) edges are depicted by dotted (dashed) lines. 
Due to the finite pulse length a, the actual spectral distribution reaches over these edges by 
0(a), where a = 7 2 (1 + f3) 2 /a. The gray shaded areas mark the overlapping regions with a 
width of 25". 

The generation of the subpeaks can be described as an interference effect [37]. Thus, when 
the harmonics overlap, for a fixed value of u' there are contributions from different harmon- 
ics and their interference reacts very sensitive to subtle changes in the phase of the jzf^ 1 
functions, see Eqs. (50) and (51). The difference in the spectral distributions looks qualita- 
tively similar to figure 1 of [38], where the influence of classical radiation reaction force on 
the spectrum was studied. The radiation reaction force also provides an electron recoil in 
the classical calculation, slightly changing the phases and leading to a modified spectrum. 
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V. SUMMARY 



In this paper, we discussed non-linear Compton scattering in the Furry picture and employed 
light cone coordinates for temporally shaped laser pulses. We emphasized the structure of 
the Volkov wavefunctions in a pulsed laser field. We clearly uncovered the modifications of 
the electron wavefunction due to their interaction with the laser field. The S matrix element 
for non-linear Compton scattering was evaluated in the framework of Volkov states within 
the Furry picture. An expression for the cross section which is independent of the pulse 
shape and pulse length was presented. 

We focused on the differences between classical calculations of non-linear Thomson scat- 
tering and quantum calculations of non-linear Compton scattering. In both cases, spectral 
broadening and harmonic substructures have been found, which are, however, shifted in the 
quantum calculation. These harmonic substructures, still lacking an experimental verifica- 
tion, are interpreted as an interference effect. We found that for ye <C 1 the differential 
quantum transition probability in many cases coincides with the classical Thomson scat- 
tering result also for pulsed laser fields, i.e. quantum effects are mostly negligible in this 
regime. Also the total Compton cross section coincides with the Thomson cross section a T 
for y e < 1(T 2 . 

As a main result, we presented a scaling law, connecting the classical and quantum spectral 
densities for arbitrary yg, e.g. yp > 1CT 2 , relating the classical and quantum results. The 
remarkable feature is, that the substructures of individual harmonics are also simply scaled. 
One might speculate that the scaling law may also be applied for arbitrary laser beams, in 
particular, for strongly focused beams. Hence, it might serve as a tool to add recoil effects 
to results obtained within classical Thomson scattering models. 

Furthermore, we also found regions in phase space, where the differential probabilities for 
Thomson and Compton are different, although yp < 10~ 2 and the total cross sections co- 
incide. This happens for sufficiently large values of a , when the individual harmonics are 
overlapping. In these regions of phase space, both spectral densities show different, almost 
erratic behavior. This observation is in qualitative agreement with previous studies of the 
effect of the radiation reaction force on the spectrum of non-linear Thomson scattering. The 
radiation reaction force introduces an electron recoil in Thomson scattering to the classi- 
cal picture. Of course, the scaling law is not applicable in the regions where the spectral 
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densities show this erratic behavior. 
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Appendix A: Light cone coordinates 



1. Basic definitions 

The light cone components of a four-vector x^ = (t,x,y,z) are defined as £±,x_i_ by 

£_ = t — Z, X + =t + Z, Xj_ = (x, y). (Al) 

In these new coordinates, the scalar product reads x • y — \ (x + y- + X-y+) — xj_ • yj_. 
Arranging these components as a four- vector = (x + , x_ , xj_ ) , one may introduce a non- 
diagonal metric 



/ | \ 

h o o o 



(A2) 



0-10 
\ -1 J 

with y/—g = 1/2, thus d 4 x — > \ dx + dx_ d 2 xj_. The inverse transformation is given by 

Xq = }^{X-\- -\- X— ), = ^(^+ — 

More generally, light cone coordinates can be introduced by projecting onto a light-like 
four- vector, which is in our case given by the laser wave vector = um+ = u>(l, n) defining 
a direction n with n • n = 1. The coordinate system will always be aligned such that 
n = (0, 0, —1). Then the light cone components of a vector may be defined by projection 
via 



B + = B-n + , B_ = B-n^, 
B ± = B - (B n)n 



(A3) 
(A4) 



with n^L = (1, — n). Additionally, we need the orthonormal transverse basis vectors e«,i = 
(1,2) with €i ■ n± = and • €j = — i.e. the transverse part of £> M may also be defined as 
Bj_ = (Bi, B 2 ), Bi = B ■ 6j. With this convention of light cone coordinates, the conjugate 
momentum to x + is P_ and vice versa. 

2. Kinematics in light cone coordinates 

Special light cone coordinates may be defined with respect to the laser momentum = 
w(l,n) with n = — e z . Let be the momentum of the incoming electron = (E p ,p) = 
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7717(1, /3) = mu^ with /3 = — f3n and the momentum of the outgoing photon k'^ = (u/, k') = 
u/(l,n') = u/n' M with n' = (cos </? sin 6>, sin (p sin 6>, cos 9). Then the light cone components of 
these vectors read 

= 2u, k + = 0, k ± = 0, 

k'_ = u/(l + cos#), £;+ = u/(l - cos#), k' ± = J sin 0(cos</?, siny?), (A5) 
p_ = 7717(1-/3), p+ = 7777(1 + /?), P± = 0. 

The free particle dispersion relation £^ = p 2 + m 2 reads in light cone coordinates 

p. = Pi+m2 . (A6) 
Using momentum conservation, one obtains for the components of p' 



(A7) 



P± = P± - k l> 
P'+ = P+ ~ K- 

It is worth noting that there is in general no conservation law for the component p'_. That 
component is fixed by the dispersion relation (A6) yielding 

m 2 + if/ J' m' + (p± - k/) 2 ,.„, 

All the components of p' have to be considered as a function of u' and n'. The projection 
of p' onto k may be rewritten as 

k ■ p' = k_p' + = k_(p + - k' + ) = k-p- k- k', (A9) 



from which the important relation 

1 1 k ■ k' 



(A10) 



k-p k-p' k ■ p k ■ p' 
may be derived. Eventually, we present a useful collection of relations between momentum 
components: 

,ii 2p • k! v • k! , . „ „ . 

k'_+p'_-p- = = 2^—, All 

p + p ■ k 

2p'-k' p + (k'_ + p'_ - p_) p + 2p-k' 
V = 2~ = 2 = ~ 2~> A12 

p' ■ k' = P± = P'k , . 

p ■ k' p' + p' ■ ' 

k ' k p ' k / a -i j \ 

X = 1 r = v^J- (A14) 
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Appendix B: Fourier expansion of the Volkov wavefunction 



Here, we provide a Fourier expansion of the Volkov state (5) in a pulsed laser field to gain 
further insight to the structure of the Volkov state and to use it as an alternative way to 
obtain the matrix element for non-linear Compton scattering as a convolution in momentum 
space. We define the Fourier transformations in X-,x + and xj_ directions with different signs 
according to the metric (A2) 



F(Q_l, Q+, Q-) = § / d 2 x ± dx_ dx+F(x ± , x-, x + )e L 2 {Q-*++Q+*-)-iQ±™± , 

The Fourier transform of the Volkov matrix C p (x) (see Eq. (11)) reads 

C P (Q ± ,Q + ,Q_) = (27r) 3 5 2 (Q ± -p ± )5(Q + -p + ) 

x [%(Q-) + <y^_^(Q_) + ^MQ-)] } 



(Bl) 



(B2) 



with 



(Q-) = J dx + e^—P-^+g\ N \(x + )e iNuJX++lfix+) . (B3) 



The functions describe the nontrivial, pulse dependent momentum distribution of the 
Volkov wavefunction. Using the Fourier representation of the Volkov state, the matrix 
element for non-linear Compton scattering (13) reads 

S fl = N (27rf5 2 (p' ± + k' ± -p ± )5(p' + + k' + - P+ )x 

1 (B4) 

with the convolution 

{% *<$ m ){k>_) ee | - Q_)Sf m (Q_). (B5) 

Comparing (B4) with (26) and (27) we find 

< = % * %, <\ = % * + ^ * Sf , 

(B6) 
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The representation (B4) of the Compton amplitude together with (B5) furnishes another 
interpretation of the subpeaks in the Compton rate as the overlap of the momentum distri- 
butions of the incoming and outgoing Volkov states. 
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